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Abstract 

We consider an ideal gas of infinitely rigid rods near a perfectly repulsive wall, 

and show that the interfacial tension of a surface with rods on one side is lower 

when the surface bends towards the rods. Surprisingly we find that rods on 

both sides of surfaces also lower the energy when the surface bends. We 

compute the partition functions of rods confined to spherical and cylindrical 

open shells, and conclude that spherical shells repel rods, whereas cylindrical 

shells (for thickness of the shell on the order of the rod-length) attract them. 

The role of flexibility is investigated by considering chains composed of two 

rigid segments. 
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I. INTRODUCTION 



Macromolecular surface interactions play an important role in a variety of phenomena: 
stabilization/destabilization of colloidal solutions [l|], where the "surface" is the surface of 
a colloidal particle; surface effects in liquid crystal displays 0, where the "surface" is the 
walls of the container; and modification of bending rigidities of bilayers ||, in which case 
the "surface" is a membrane. 

Both rigid and flexible macromolecules can induce a modification of the properties of 
interfaces, either by adsorption or depletion; here we concentrate on the surface physics 
of non- adsorbing, rigid and rod-shaped macromolecules. Interest in solutions of rigid rod- 
like objects can be traced to the work of Stanley who first extracted the tobacco mosaic 
virus (TMV) @] and to the subsequent observation of the nematic phase in TMV solu- 
tions ||, which was later explained by the seminal theory of Onsager ||. Many other 
rod-like molecules exist in the biological realm, ranging from DNA in its a-helix form to 
fibrils of amyloid /3-protein, the molecular agent at the origin of the Alzheimer disease ||. 
Solutions of rod-like molecules in a large variety of mineral and organic systems have also 
been studied [01. 



In a pioneering study of surface interactions of macromolecules Asakura and Oosawa |T0 
showed that the steric depletion of rigid macromolecules at a flat surface increases the 
interfacial energy of the system, implying that two surfaces in a solution of non-interacting 
rods will attract one another when they are at separations smaller than the rod-length. 
This effect, which is at the origin of colloidal destabilization, is also present in the steric 
depletion of spherical particles, or particles with other shapes - but with a different range 
and magnitude. Asakura and Oosawa theory has since been extended to include effects of 



nematic ordering in the bulk [11] or effects of rod- rod excluded volume interactions [12 



in the latter work the authors show that when the modification of the rod-rod excluded 
volume interactions close to the surface is taken into account, a rod solution may stabilize a 
colloidal suspension, rather than destabilizing it. However, almost all the available studies 
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on rod-surface interactions are devoted to flat geometries, except for Auvray's work fl3 
where the author makes a first attempt at considering curved geometries; we use some of 
his constructions in our work, which extends the original Asakura and Oosawa theory to 
include effects of surface curvature in the contribution of the rod depletion layer to the 
surface tension. For flexible surfaces, e.g. membranes, in rod solutions, our results allow 
us to derive the renormalization of the bending rigidities of the membrane. Some of our 
findings are qualitatively different from the corresponding results for flexible polymers in 
solution. This prompted us to consider two rigid rod-segments joined in the middle as a 
first step in studying the cross-over from rigid rods to flexible polymers - rigid rods being 
a polymer of which the persistence length is longer than its length, and the rod-segments 
joined in the middle being one with a persistence length equal to half of its length. Not 
surprisingly we find that when one allows for a hinge in the rod, the results are closer to the 
flexible polymer results. 

Even though steric effects are purely entropic and exist for any non-adsorbing macro- 
molecular solution, the strength of the effects is larger for rods with large aspect ratios. The 
following simple comparison of the order of magnitudes of the contributions to the surface 
tension with those caused by spherical particles illustrates this. The typical scale of the 
energy density in a solution of colloidal spheres of radius tq and particle number density p& 
is pb- Integration over the layer thickness results in corrections to the interfacial energy 
of the order of A7 ~ ksTp^Q. On the other hand, the scale for the interfacial tension in 
most liquids is of the order of 70 — ksT/a 2 where a is a microscopic size. For instance for 
a ~ O.lnm, 70 is of the order of tens of mN/m. The depletion of spherical particles, even 
at order unity volume fractions = pb(47r/3)rQ, induces thus corrections to the interfacial 
tension which are a factor (a/ro) 2 lower than typical values. Curvature corrections to the 
interfacial energy are of the form: A7 ~ ksTp^r^l + Cir /R + C^r^/i? 2 ), where R is the 
radius of curvature and C\, C2 are two numerical constants. This sets the order of magnitude 
of modifications of the bare elastic constants: An ~ ksTp^rQ. Even at the upper concen- 



tration limit pb ~ l/ r o> these corrections are at the lower end of the range [14l, 1 — 20/csT, 



of most bare elastic constants. For an isotropic rod solution the upper concentration limit 
is the Onsager concentration p\ = 4.2/ (L 2 D) ||, where L is the length of the rod and D its 
diameter. Now interfacial tension contributions are of the order A7 ~ ksTpbL, but even for 
rod diameters of order of the microscopic length a, and at the Onsager concentration this 
contribution is still a factor (a/L) smaller than typical interfacial tension values. However, 
modifications of the elastic constants are here of order An ~ ksTpbL 3 . At the Onsager con- 
centration this is a factor (L/D) larger than ksT, and therefore, may lead to a substantial 
modification of the elastic constants of even rather rigid phospholipid membranes with elas- 
tic constants as large as 20 ksT, even at low rod concentrations where rod-rod interactions 
are negligible. 

In the next section we introduce the basic thermodynamics necessary for our discussion; 



sections |T| and [IV| contain our results for rods close to surfaces, and rods in shells, respec- 
tively. We make an initial attempt at understanding the cross-over from flexible polymers 
to rigid rods in section [V]; the last section is devoted to a discussion and summary of our 
main results. 

II. THERMODYNAMICS 

We consider an ideal gas of rods of length L in the presence of surfaces which repel 
the rods, in several different geometries. The thickness of the rods is taken to be zero, i.e. 
much smaller than all other length scales in the problem. We parametrize the possible rod 
configurations by the coordinate f, and two angles specifying in which direction the rod 
points, uj = uj(6,(p). Kinetics is ignored throughout the treatment^ 

The relevant potential describing the thermodynamics of the system can be written as, 
(k B T = 1): 

Q[p(f,uj)]= / dr / dujp(f,Co) [log(v p(f,u))/e) - (p b - U ext (r,u))} (1) 



1 Kinetic effects will be smaller by a factor of (thermal wavelength/L) ~ 0. 
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where v is some normalization volume, pb is the solution chemical potential, and U ex t is 
the hard wall interaction potential that is either infinite or zero, depending on whether 
the configuration of the rod is allowed by the 'hard wall' requirement or not. Functional 
minimization of Q with respect to p(f, ul) gives the equilibrium density profile: 

v Air 

with pj, the rod number density of the bulk solution. The local, position dependent, number 
density of particles p(f) is simply the sum over all allowed angular configurations of p(f, ul). 
The excess surface energy is then calculated as follows: 

A ^ = n[p(f^)]-n[ Pb /(An)] =p fdr r dw^ _ ^ UesM) (3) 
S J S J Aix 

where S is the surface area and the volume integral runs over the space available to the 
solution. 



III. ONE SURFACE IN A ROD SOLUTION 

Here we consider surfaces with rods on one side only. Results for the case with rods on 
both sides can be easily found by adding the results for the inside and the outside. Our 
calculations ignore fringe effects that are unimportant in the thermodynamic limit, e.g. the 
flat wall is treated as translationally invariant. 

We first rewrite equation ([|) in pressure units (ksTpb = 1) and explicitly state its rod- 
length dependence 

A 7 (£)=/f J§ h(r,u,L) (4) 

If we represent the space coordinates of the rod by the position of one extremity (see fig- 
ure ([]])), the function h(f,uj,L) is unity when a rod of length L, starting at r and point- 
ing in the direction uj intersects the surface; it is otherwise. This can be written as 
h(f, uj, L) = H(L—(s—f)-u}), with H{x) the usual step function H(x > 0) = 1 ; H{x < 0) = 
and s is defined as follows: The half-ray orginating at r and travelling in the direction of u 
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may or may not intersect the surface. If it does not then we do not have to define s, since 
all such configurations give H = 0, and thus do not affect any of the below expressions. If 
the half-ray does intersect the surface, we define s as the intersection point that is closest 
to r. Notice that s is independent of L. 

It is convenient to calculate the L-derivative of the excess surface energy: 

with 6g the angle between the rod and the normal to the surface at point s and g(s, u, L) 
a function restraining the angular integration space for a rod with one extremity at s. The 
integral over (is ranges over the whole surface. The factor cos 6^ is the appropriate Jacobian. 



A. Convex bodies 

For any surface delimiting a convex body, the function g(s,u,L) restricts the angular 
integral to half of the angular space, independently of the actual surface shape. 

Q I r2n «r/2 1 

— A 7 (L) = — / d(j) / dOsine cos9 = - (6) 
oL 47r Jo Jo 4 

Integrating equation (g) with the boundary condition A^(L = 0) = one gets the three 
dimensional result for the surface tension of convex bodies immersed in a rod solution: 

A T (£) = \ (7) 
This result can easily be extended to convex bodies in an arbitrary dimension T>: 

A 7 (L) = c(V) L (8) 

with 

rg/2) (P-2)H J \ ; for V odd 

1 1 2 v /?r( CD + D/2) (D _ 1)n I i . for v even < J) 

In particular, c(2) = 1/rr and c(3) = 1/4. 



B. General Surfaces 



For a general surface, the function g(s,u,L) might restrict the angular integration in 
equation ([5]) to less than half of the angular space. In this case we can still integrate over 
half of the angular space and further subtract the non-available angular sector. In three 
dimensions: 

A A7 (L) = \-J§ / ^ cos9,H{L - (P - s).Q) (10) 

where s' is another surface vector. Differentiating with respect to L and performing the 
angular integration leads to 

d 2 A ,-rs f ds f dl n „ . N 

A 7 (L) = - — —— cos ^ cot % (11) 



dL 2 " ' J S J AirL 2 
The vector i defines all the points on the surface at a distance L from a given point s. In this 

three dimensional case this is usually a one dimensional curve. The angle 8g is measured 

with respect to the normal of the surface at point s and the angle Of is measured with 

respect to the surface normal at point t. Contrary to the convex case, the excess surface 

energy of a non-convex body immersed in a rod solution does depend on the actual shape of 



the body. Moreover, because the contribution of the second order derivative ( HJ) is always 
negative, any long- wavelength concavity will reduce the surface tension of a surface with rods 
on one side or both sides. Therefore entropic effects in such systems will render flat surfaces 
unstable towards long-wavelength curvature! It also implies that for a fixed geometry the 
surface tension is a concave function of the rod-length, L, and will have the general form 
depicted in figure (§). 

Equation (|ll|) can easily be extended to spaces of different dimensions. For instance, in 
two dimensions, where the "surface" is actually a one dimensional curve of perimeter P, one 
has 

d 2 



dL 



A 7 (L) = -Jp E^ cos d i™t e 9 (12) 

The sum runs over all (usually two) points q at a distance L from a given point I on the 
surface. 



C. Some examples 



With the above methods, Aj(L) can be calculated exactly for the inside and outside of 
a circle, a sphere, and a cylinder, each with radius R; and also for two flat walls separated 
by 5, i.e. a gap. We would like to emphasize that all of these results can also be derived in 
an independent way which is also exact, namely by "counting beans", i.e. by just doing the 
integrals in Eq. (|]) straightforwardly (see section [IIIE| below). We have explicitly checked 
that all of the results quoted here agree with these straightforward calculations - especially 
in the cylinder's case both calculations are quite involved algebraically, and this check is 
absolutely non-trivial. Of course, the method of computing the second derivative and then 
integrating up is the shortest route to these answers. 

We now define e = L/R; this is for convenience only and does not mean that we take e 
to be small unless specifically stated. 

We first calculate results for two dimensions. Inside of a circle we have, (see figure (Of)): 



cos0£ = e/2, cot 6^ = cot#£ = (e/2)/ ^1 — e 2 /4^, and the sum over q amounts to a factor 
of two. Thus: 

|L A7 » = _^ cosW = _£_^ (13) 

Using A 7 (L = 0) = 0, and ^ A 7 (L = 0) = ± we find: 

Ai in 1 n 1 / 

-JfL = - sin" 1 (e/2) + ±-{\ - eV4 (14) 
For three dimensions we consider rods inside spheres and cylinders. Inside of a sphere we 



have, (see figure (|)): cos6> s = e/2, cot 9 e = cot6 s = (e/2)/ (Jl- e 2 /4j and I = 2nL sin6> s . 
Thus: 

J^ A 7i = ~ *™0s cos0 s cot9 s = -il (15) 



It follows that: 



_ 1 1 2 



4 48 f < 16 > 



The calculation for the inside of a cylinder is a bit tedious, but also doable. We find: 

d 2 



dL 2 A ^cyl 3nRe3 



in 



e 2 " 



e 2 



(17) 



where the functions E and K are the complete elliptic integral of the first and second kind 
respectively. The cylinder result can in principle be integrated to give the special functions 
MeijerG. The straight-forward "bean-counting" mentioned above leads also to an integral 
form that can only be done explicitly as a perturbation series in e. We checked that the 
coefficients in this series agree with the expression for the second derivative in Eq. (|17])- The 
perturbative result for the inside of the cylinder is: 



M%1 _ 1 1 2 



4 m 



e 2 + ■ ■ ■ (18) 



The case of the flat gap has a different topology, and therefore is worth doing as an 
example, (see figure fl5|).) For L < 5, clearly there are no configurations of a rod that bridge 
two points on the surface, and therefore: A7 gap = L/2, (to assure consistency with sections 
to come we took S to be one side's area.) Whereas for L > 5, we have: 9i = 6 2 = 6, 
cos 8 = S/L, and £ = 2nL sin#. Thus: 

d 2 . cos 2 6 5 2 

A 7gap = r~ = -7J (19) 



It follows that: 



8L 2 ,9ap L L 3 



Algap J 5/L-5 2 /(2L 2 ) ; L>5 



(20) 

1/2 ; L < 5 

This is to be compared to Eq. (|25| ) below. 

We also developed a method to compute the small L expansion of the surface tension for 
a general surface in two and three dimensions. We leave this discussion to appendix and 
just note here that the expansions of the expressions for the circle, sphere, and the cylinder 
agree with this general method's predictions. In particular, that A7 is an odd function of 
L for sufficiently small L follows from this expansion. 
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D. Random Insertion Sampling 



We also tested the surprising prediction of equation^ in the case of an ellipsoid, the 
simplest non-trivial case in 3 dimensions, by performing random insertion sampling. Since 
the prediction for the surface tension is in the limit of a non-interaction gas of rods, we need 
only consider the behavior of one rod of length L in our numerical experiment. 

The algorithm is as follows. We choose a volume Vo around the ellipsoid of semi-axis 
A, B, C that is large enough to encompass the ellipsoid and a shell of thickness L/2 around 
the ellipsoid. We then randomly select M center of mass positions and orientations for the 
rod of length L. The center of mass must be in Vq excluding V e , the volume of the ellipsoid 
itself. By noting the fraction / of insertions that lead to overlap between the rod and the 
ellipsoid we can express the surface-tension as: 

A 7 (£) = / • (21) 

in units of ksT 'p b - S e denotes the surface area of the ellipsoid. 

When we choose M = 10 8 , our numerical experiment confirms the prediction Aj(L)/L = 
1/4 within an error-bar of 0.5 * 10~ 4 for a wide variety of choices for A, B, C, and L. The 
error-bar on the fraction / is J ^rp^ 



E. Comments 

Although the method of "bean-counting" is less efficient than the one described above, it 
is nevertheless useful for gaining more insight into the reasons for the counterintuitive results 
above. In this method one performs first the angular integration over p(f, u) to obtain the 
rod density p(r). In the particular cases of flat, cylindrical and spherical surfaces, the excess 
surface energy is then obtained from 

A 7 = | dz(p b -p(z))J(z,R) (22) 

with z the perpendicular distance from the surface, and J(z, R) the Jacobian, which depends 
on the geometry. For flat surfaces J(z,R) = 1; for the outside of a cylinder of radius R, 
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J(z,R) = 1 + z/R, and for a sphere J(z,R) = (1 + z/R) 2 . For these three geometries 
there is only a one-dimensional integral (over z) to be carried out, once p(z) is determined. 
Differently from the previous sections, we define hereafter the rod coordinates by the position 
of its center of mass. 

For example, for a plane with rods on one side (see figure (^j)), we have: 



rL/2 

A 7/ = / dz 







I r2n rn— arccos(2z/L) 

1 - — / d(j) / sin e d6 

4 7T JO J a.rccos(2z / L) 



rL/2 


r 2zi 


= dz 




Jo 





4 ^ 



Notice that the rod concentration p = p^lzjL is different from the monomer concentration. 
The evaluation of this second quantity is slightly more involved but can be calculated fol- 
lowing Auvray lfL3| . This gives: Cf(z) = Cb (z/L)(l —log(z/L)), where Cf, = pbN is the bulk 



monomer concentration for a solution of rods with N monomers. 

The calculation of the phase space available to a rod in the presence of a curved surface 
is somewhat different from the one for a flat surface, for now it is possible for the rod to 
contact the wall in two different ways (see figure ([?[)): when the rod is sufficiently far from 
the surface the tip touches the wall upon rotation, whereas when the rod is sufficiently close 
to the wall, it will be tangent to the surface at some point along its length. The critical z 
value separating these two regions is z™* ph = yj (L/2) 2 + R 2 — R. The phase spaces arising 
from these two configurations have different functional forms. In the inside of the sphere, 
of course, the tangent construction is not possible, but now there is a minimum distance 
beyond which it is not possible to place the center of mass of the rod - this happens at 
z *!sph = \] R 2 ~ (L/2) 2 — R. A closer look at the results from the previous sections reveals 
that the surface tension is non-analytic in the curvature; e.g. that even for arbitrarily 
small e = L/R, one has neither A^ l h (e) = A^ h (-e) nor A^ h (e) = -Ay£ h (-e). This 
invalidates the usual extrapolation of the free energy results from the cylindrical and the 
spherical shapes to arbitrary bending states. We can see that the difference between z^f ph , 
and z™ sph is the basic reason for this non-analyticity. It is also possible to compute the 
monomer densities in the spherical configuration, we present these in appendix [FJ. 

From equation ( p3l) it is clear that differences between the excess energy of a flat and a 
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curved surface depend on two factors. The first is the configurational part measured by the 
differences in the rod density profiles, the second one is associated with the space available to 
the center of mass, and is measured by J(z, R). In the case of hard sphere solutions, where 
there is no coupling between configuration and curvature, only J(z, R) is responsible for 
energy differences; the profile of particle concentration, p(z) = pb Q(z — r ), is independent 
of geometry. It is easy to show that in that case one has A7 = /cbT/vo(1 + r /(2i?)) for 
cylindrical surfaces and A7 = fcsTpfer (l + r /R + r\j (3i? 2 )) for spherical ones. In the case 
of rod solutions, there is a coupling between the configuration of the rods and the surface 
curvature: the concentration profiles are a function of the surface curvature. Given this 
coupling, the exact correspondence among the free energies of all convex-shaped rod-systems, 
proven in section |III A| , is rather surprising. This implies that curvature contributions to free 
energy fl22|) from both the profile and the Jacobian integrate out to a total vanishing value. 
This cancellation is non-local, i.e. at a given z away from the surface the total phase space 
for the center of mass does not change as to compensate for the configurational changes 
in the profile. One can gain further insight into the issue by recalling equation (p2"D, and 
considering figure (^j). The surface tension is just the total "area" between the curves p(z) 
and pb weighted by the appropriate Jacobian. The surface tension for the plane is simply the 
sum of the shaded areas I and II, whereas the surface tension for the outside of the sphere 
is the weighted "area" of I only, which happens to be the same. The relevant weighted areas 
are shown in figure ([|). 

Also, we pointed out that in general the surface tension decreases as the surface develops 
concave regions, e.g. ^Ysph < ^7/> i- e - the sum of the weighted "areas" of the regions I, II, 
and III in figure (|8]) is smaller than the area of I. Therefore a flat surface with rods on one 
side would be able to lower its energy by bending towards the rods in a spherical shape; and 
if there are rods on both sides, since Aj^h < A7J, they would destabilize the flat surface 
towards bending. 
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F. Bending Rigidities 



For a flexible surface, e.g. a membrane, one computes the renormalization of the bend- 
ing rigidities by assuming that A/ c , the excess free energy per unit area of the system is 
expandable in a power series in the curvature, i.e. the excess free energy is analytic in the 
curvature. Then up to quadratic order in the principal curvatures and j^, this expansion 
can be written in terms of the mean and Gaussian curvatures of the surface [0]. In terms of 
the principal radii, the mean curvature is H = |(^- + and the Gaussian curvature is 
K = A general form of Af c to this order is given by: f c = 2Ak(H — c ) 2 + AkK, 

where Ak and Ak are the contributions from the excess free energy to k and K, the bare 
elastic constants of the membrane, cq is the induced spontaneous curvature of the system. 

As an example we calculate first the modifications of the elastic constants of a membrane 
immersed in a dilute solution of hard spheres. As quoted above A7 = ksTpbro^l + r / (2R)) 
for cylindrical surfaces (l/i^ = 1/R and 1/R 2 = 0) and A7 = k B Tp b r (l +r /R + r%/{3R 2 )) 
for spherical surfaces (1/Ri = I/R2 = ^/R)- This leads to Ak = 0, Ak = (2/3) k B T p b r\ 
and Co = 0. 

For a membrane exposed to a rod solution, the non-analyticity of the free energy might 
imply that the extracted bending rigidities are not meaningful for all states of curvature 
of the membrane. In particular, bending states where the principal radii of curvature have 
opposite signs might not be described by the usual Helfrich expansion. With this note of 
caution in mind, we computed the renormalization of the bending rigidities for a membrane 
exposed to a rod solution on both sides: 



c = 

(24) 



Ak = -i ( P6 L3) = §^1; 



Ak = +^ (p b L 3 ) = J^t5 



It is interesting that there is no spontaneous curvature induced by the rods, since there 
is no symmetry reason why this should be the case. Though, our "small L analysis" in 
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appendix [A] shows that this is an (almost) general property of rod-systems. The modifica- 
tions of the elastic constants (a decrease of k and an increase of K indicating a preference to 
form periodic minimal surfaces) are similar to results for depletion of other macromolecular 
species. However, the magnitude of this effect can be here much larger than fc^T, a feature 
not found in other depletion systems. It is of order of (pbL 3 ) ~ (Pb/pt) (L/D), where D 
is the thickness of the rod and p* b is the Onsager concentration for rigid rods. Even at the 
Onsager concentration the enhancement factor L/D can lead to very large contributions. 



IV. TWO SURFACES IN A ROD SOLUTION 

Here we consider rods in shells where both surfaces of the shells are repulsive. The shells 
are taken to be in equilibrium with a bulk reservoir of a rod solution at bulk concentration 
Pb- The thickness of the shells is denoted by 5, the perpendicular distance from the inner 
surface by z. The volume, and the thickness of the box are kept fixed as it gets bent. 
When 5 > L, the results in this section follow from the ones in the previous sections, since 
the effects of the two surfaces do not interfere with one another. When 5 < L however, 
there is interference and one needs to calculate the allowed phase space by accounting for 
interactions with both surfaces. 

We consider rods confined inside flat or curved shells, but in equilibrium with a rod 
solution. By comparing the energies of flat, spherical and cylindrical shells, we hope to 
investigate for rods a question that has recently been investigated for flexible polymers fl5| 



and mesoscopic quantum systems [jl6[| : do the curved regions in confined systems act as 
attractive potentials, where the concentration of particles will be greater, or as depletive 
potentials from where particles will be repelled ? Table p] summarizes the systems studied 
so far, including the case of confined rods treated below. 
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A. Flat Strip 



A straightforward calculation for a flat strip of thickness 5 = y L confining a rods solution 
yields: 

A 7 f \y-\y 2 ]y <l 

> (25) 

\ ; y> 1 

As usual, the interfacial energy is expressed in units of k B T p b . Note that for separations 
larger than the rod-length the interfacial energy is just the sum of two contributions from 
the inside surfaces. When the surfaces are taken to be also in contact with a solution outside 
(see figure (.|10J)) we get: 

A 7 J^_ J \ + v-\v 2 ; y<i 

y > i 



B. Spherical Shell 

We are able to calculate the available phase space, and thus the surface tension of a 
spherical shell exactly. The result has four regions, across the boundaries of which the 
functional form of A7 changes. This happens because one needs to consider the relative 
positions of z^f, zl n s , which were defined earlier, and z , the point separating the two 
regions where the rod only touches the upper surface and the one where it only touches the 
lower one. The extent of the first two regions vanishes with curvature; the results below 
tend asymptotically to the ones for the flat strip in the limit of small curvature. We took 
S as the surface area of the flat box with the same volume and thickness. We leave the full 
results to the appendix and give here the expansions in e = L/R. 
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A 7 , 



m 

sph.sh. 



(i-¥ 



48 y 



+ 0(e) f 



l + ^ 2 y 2 -i) + o(e 



y<^f¥^- 1 7 



i 4- i 

4 "I" e 2 



< i/ < v 1 
l+£-I<y<l 

y > i 



(27) 



One can check by comparing Eq.(|25D with Eq.(|27|) that for any separation the surface 
tension of the spherical shell is larger than that of the flat box. We conclude therefore that 
rods prefer to stay in flat regions of the confining box rather than in spherically curved ones. 
This can be further checked by computing the average rod concentration inside the flat shell 
and comparing it to the value inside the spherically bent shell; the latter will be smaller than 
the former. We show below that cylindrical shells induce the opposite, i.e. attract rods, for 
L less than but close to 5. Results for a spherical shell immersed in a a rod solution can be 



easily obtained by adding to Eq.(p7|) the outside and inside contributions of Eq.fllED 



C. Cylindrical Shell 



The cylindrical shell calculation is substantially more involved than the ones for planar 
or spherical shells. We resorted to a perturbation in (1/R) in this section, and stayed at 
separations close to the rod-length, i.e. the results below for 5 < L should be understood to 
be valid only in a finite region close to L. Again, choosing S as the flat box's surface area, 
we find: 



A^H 1 -^^ 1 - 4 ^ 2 ^'--^ 1 (28) 

I _ J_ e 2 , . . . ■ V>1 

The results for the cylindrical shells are, at least for separations close to the rod-length, 
opposite of those for the spherical shells: the surface energy in this case is smaller than that 
for the flat strip. Rod-like molecules will therefore be attracted by cylindrically bent shells. 



16 



D. Bending Rigidities 



As explained above, for small curvatures, it is possible to describe the elastic properties 
of a surface, by three parameters: the spontaneous curvature cq, and the bending rigidities 
k, and R. In this paragraph we compute these parameters for a flexible box, such as a 
membrane, exposed to rods. For all cases Cq = 0, as it should be due to the symmetry of 
the geometry. 

For an open shell with rods on the inside only we find: 

5 V (29) 

1 (_1 + l 2y 2 _ 2y 4) . y<1 



pi d 



22.9 V 1 » "*>> ' » ^ * 

( ^9 (1 + % 2 ) ; y > 1 



For a strip immersed in a solution of rods, i.e. a system with rods inside and outside the 
shell we find: 



PI 



Ak = _^l jsH + V-V);i/<i 

h ]V>1 (30) 

r± l \ h (- 1 + 20 ?/ 2 - V) ; y<i 



(1 + 16?/ 2 ) ; y > 1 



22.9 



V. FROM RIGID RODS TO FLEXIBLE POLYMERS 

A rigid rod can be viewed as the particular limit of a freely-hinged polymer chain with 
one monomer only, see figure (|TTD. In the reverse limit, when the polymerization number N 
is very large, the equilibrium configurations of such a chain are well described by Gaussian 
statistics. We now describe to some extent the crossover between these two limits. 

For a membrane in contact with a dilute solution of Gaussian polymers, it is possible to 
calculate the excess surface energy [17]. For instance, a sphere in T> dimensions immersed 



in a such a solution increases its surface energy by 
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2 ATl/2 , D-liVa 2 , (D- !)(£>- 3) iV 3 / 2 



3 



A7^(iV) = 4= « + ^ ^ + 1 "T" ^ -^T- + ■ ■ (31) 



2 R 6v^F i? 2 

with a the monomer size (and total length L = Na). For a sphere in three dimensions this 



reduces to: 



A-fs P h, P oiy(N) = -|= N 1 ' 2 a + ^ + • - • (32) 



Note that for the Gaussian polymer case, there is an analytical expansion of the free-energy 
in powers of the curvature, the case of a solution inside a sphere can simply be obtained 
from the result above by a sign permutation 1/R —1/R. The scaling of the first term in 
equation fl3"2"|) implies that for a plane the excess surface energy decreases by a factor of y/2 
when one increases by a factor 2 the number of links in a chain of length L. 

In the rod limit (N = 1), one has in the flat case, A7/(JV = 1) = (1/4)L. The (N = 2)- 
calculation for the flat case is straightforward and gives: 

A V (N = 2) = (33) 

where L = 2Na is still the total length of the polymer. Interestingly, even for such a small 
monomer number, doubling the number of monomers reduces the surface tension by 1.67, a 
factor not too different from \/2 in the Gaussian (large N) limit. It is also possible to write 
down a general integral expression for the surface tension near a flat wall for any N where 
all the integrals are in principle doable; we left this to the appendix, Eq. (|CT|). 

The surface tension for (N = 2) in- and outside of a sphere can be calculated using the 
same methods as for the rods (N = 1). We leave the exact outside result to the appendix, 
Eq. (|C2|) , and expand it here in the curvature: 

A 7 ^(iV = 2) = — + — e- — e 2 + --- (34) 

IsphK ) 24 g 4 24 g V I 

The exact inside result is simply: 

A 7 ™ (N = 2) = — - — e - — e 2 (35) 

lsph\ ) 24 g 4 lg20 V ) 
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Notice that there is now a spontaneous curvature. Moreover, the sign of it is positive for the 
outside, as it is for ideal polymers ||17|| . Also, the asymmetry between the in- and the outside 
cases, which is the cause of the non-analyticity in the curvature free energy, is exactly one- 
third of what it is for (N = 1). Furthermore, the overall magnitude of the e 2 term is about 
half of what it is for the rigid rod. This term does not exist for the polymer. Therefore, 
even the introduction of one single link in the rigid rod, brings the results closer to those of 
the flexible polymer case. 

VI. CONCLUSION 

We studied the depletion interactions of curved surfaces exposed to dilute rod solutions. 
We find that in general, the excess surface energy caused by the depletion interactions 
depends on the exact shape of the surface. For instance, the surface tension of a plane 
exposed to a dilute solution of rods of length L increases by A7^ at = /csTp b L/4, with p b the 
rod number concentration of the solution, whereas the equivalent quantity for a spherical 
surface in contact with such a solution inside the sphere is AY S ph — ksTpbL/A x (1 — 
(1/12) L 2 /R 2 ), with R the radius of curvature of the sphere. 

Surprisingly, we found also that a spherical surface in contact with a rod solution outside 
of the sphere has the same excess surface energy as the flat surface A7°^ = ksTpbL/A. 
Moreover we proved that this equality holds for any convex surface exposed to the solution. 
This implies that the excess energy of any convex body immersed in a dilute rod solution 
only depends on the total surface of the body and not on its actual shape (as long as it 
is convex). A similar equality also holds in any dimension of space. For instance in two 
dimensions, the line tension for a straight line is ^fiat = kBTubL/ix with a b the "bulk" 
surface rod density, and has the same value for rods outside of a circle. 

The first practical implication of our results is that a flexible membrane immersed in a 
solution (thus exposed to the depletion rod layer on both sides) will spontaneously break 
its curvature symmetry. Also, a membrane exposed to a rod solution (on one side) will 
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spontaneously bend towards the solution. The experimental observation of these effects 
will depend on the magnitude of the symmetry breaking field. This can be measured by 
comparing the relative importance of induced modifications in the elastic constants of the 
membrane. Typical values of the bare elastic constants are in the range 1 — 20ksT. The 
induced shift to these bare constants are for most systems like hard spheres or flexible 
polymers, only a fraction of ksT. However for rod solutions, at the Onsager concentration, 
the induced modification is a factor L/D larger than ksT, and can potentially be very large. 

When the rods are confined in curved shells, we showed that regions with cylindrical 
curvature attract the rods, whereas regions with spherical curvature deplete the rods. This 
enlightens a recent discussion on confinement effects in solutions of flexible polymers fl5 



and in mesoscopic quantum systems |16|] . Table | summarizes results in a variety of systems. 

The methods that we developed for this work can be applied to other geometries and 
other situations. Work is under progress to exactly evaluate rod-mediated sphere-sphere in- 
teractions. Also of interest are the effects of attractive interactions on the bending constants 
of membranes exposed to rod solutions. 

We remark also that the steric surface interactions are identical for rigid rods and rigid 
disks near surfaces, if one takes the rod's length and the disk's diameter to be the same; 
however the disk system does not enjoy the above-mentioned L/D enhancement. The shell 
calculations also generalize to the disks if the radius of the disks is sufficiently small. 
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APPENDIX A: SMALL L EXPANSION 



When L is small (or the radii of curvature large) we can expand the excess free energy 
in powers of L. We do this first for the 2-D case. We start with Eq. ([12]). For L small, 
Yjq cos <fie cot <fi q is a function only of the point £ and the expansion of the shape around it. 
We assume that L is small enough that rods of length L starting at £ can only reach out to 
points which are nearby along the surface. If the surface is convex in this region, there are 
no points q, and the function is 0. We ignore inflection points, because they will generally 
be of measure (except for plane regions, where the function is also 0). We assume that 
there are no corners. This leaves only regions where the surface is concave. Around the 
point £ we can draw the tangent to the surface as the x-axis, with £ at x = 0, and let the 
curve's shape be given by the function f(x), (See figure (0).) The function / will have a 
Taylor expansion: 

f(x) = c 2 x 2 + c 3 x 3 + c 4 x 4 H (Al) 

where the are functions of £. We want to find the two values of x such that rods of length 



L start at the origin and hit the curve (See figure flT2|).) That is, we need solutions to 
x 2 + f(x) 2 = L 2 . There are two solutions, which for L small are near ±L. We can expand 
in powers of L, and to 4 th order find: 

x± = ±LtIc 2 L 3 -c 2 c 3 L 4 + 0(L 5 ) (A2) 

It is not hard to show that the two solutions x± are given by the same power series, except 
with the sign of the odd powers in L reversed. Another way to say this is that the two 
values of x are given by one power series in ±L. Given the shape f(x) and values x±, we 



can calculate the cosines and cotangents (see figure ([12])): 

/ xf — f 

cos (fie cot (fi q = 



x=x± 



Ll^i + wy/ \ £ 2 (i + (f 

Plugging in the expansion of x± to order L 3 we get: 



W'-fl 2 (A3) 
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cos cot (f)q = c\L 2 ±3c 3 L 3 + 0(L 4 ) 

(A4) 

E g cos0^cot0g = 2c 2 .L 2 + 0(L 4 ) 
In the sum over x±, the terms of order L 3 canceled out. Because x± is a function of ±L, 
functions of x± are also functions of ±L. The expression for cos cot (j) q above uses only 
functions of x± and L 2 = (±L) 2 , so is again a function of ±L. It follows that when we 
sum over q, all terms odd in L will cancel out in X^cos^cot^, so and A7 are odd 
functions of L for L small! If we put everything together, and integrate, we get: 

A'y(L) = an odd function of L 

(A5) 

= L ,-&M) + o{ L s) 

Here R is the local radius of curvature, and (j^) indicates that we are averaging over 
the perimeter of the object, weighting by arclength. We set jp to where the object is 
convex. The assumption that there are no corners is necessary - for example, for rods inside 
a rectangle there will be terms of order L 2 . And if L is too large, we can no longer do this 
small L expansion and A7 does not have to be an odd function of L - for example, if we are 
inside a closed body, and L is larger than all lengths of the body, A7 is simply a nonzero 
constant independent of L. 

The expansion in 3 dimensions is a little more complicated, but essentially the same. We 
again assume that there are no corners, and that inflection points outside of plane regions 
are a set of measure 0. We have a tangent xy-pl&ne and a function f(x,y). We expand 
f(x,y) in powers of x and y, starting at x 2 and y 2 and rotating out the xy term. The 
solutions to x 2 + y 2 + f(x, y) 2 = L 2 are a near-circle. Letting be the angle in the xy plane, 
and putting u 2 = x 2 + y 2 we have 

f(u,(f>) = c 2 ,ct>u 2 + c 3:4) u 3 + c 4t4> u 4 H 

C2,<f> = c 20 cos 2 + c 22 sin 2 (A6) 
c 3; ^ = c 30 cos 3 + c 3 i cos 2 sin + c 32 cos sin 2 + c 33 sin 3 
Note that c n ^ = (— ) n c n ^ +7r . We now need to look for solutions of u 2 + f(u, 0) 2 = L 2 for 
each particular 0. But this is just the 2-D problem, which has already been done. The 
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only additional complication is the arclength factor, which can be rewritten in spherical 
coordinates as: 



di = J(Ld6) 2 + (L sm9 d(f)) 2 



= #L^l-cos^ + ^(^) (A7) 
= ^L^l-(f) 8 + ^,(i(f)) 9 
The values of / at angles 4> and </> + 7r have the same power series in ±L, and L only appears 
independently in L 2 = (±L) 2 . So as in the 2-D case, A7 is an odd function of L for L small. 
The low order expansion in L gives: 

/ di cos 6 X cot 6 2 = L d<j> [c% L 2 + c 30 L 3 + 0(L 4 )] 

= L 3 f (3 c 2 + 2 c 20 c 22 + 3 c 2 2 ) + 0(L 5 ) 
Integrating twice and rewriting in terms of the principal radii of curvature, 



Once again, when taking this average it is understood that we take the quantity inside the 
parentheses to be zero if the surface is concave at the point at hand. Now in terms of the 
Gaussian and the mean curvatures: 

Arf(L) = an odd function of L 

(A10) 



= f-lf-#J^ + ^ s ) 
APPENDIX B: SPHERE AND SPHERICAL SHELL 

The probability of finding a rod's center of mass z away from a sphere can be calculated 
by first computing the allowed angular phase in figure (|7|) and then integrating over the 
angles. For convenience we define r = R/L — 1/e, and still use x = z/L. We find for these 
probabilities: 



pf P h( x )/pb 



5 Jj *,sph 



I \Jlrx + x 2 1 (r + x) 

2x + (1/4 - x 2 )/(r + x) ; x°% h < x < 1/2 

; < x < x% vh 

2x - (1/4 - x 2 )/{r - x) ■ xl n sph < x < 1/2 
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(Bl) 



The monomer concentration in- and out-side of a sphere can be computed following the 
construction in flI5| . We define A*(r, x) = 1/2 — JJr — x) 2 — r 2 + (1/2) 2 : 
For the inside of the sphere the concentration profile is given byQ: 



C sph( x ^ 



cf phJ (x) ; < x < r - ^ (1/2) 2 + r 2 



c sph,II\ x ) 



r - y (1/2) 2 + r 2 < a; < 1 
x > 1 



where 



x — x log (x) 



1 



2(r-x) \2 



|(1 -x 2 ) +x 2 log (x)) + 



(rifc) " 55^ (K A * - (1 - A *) 2 ) - (A)) 



/7 (x) = c b 



X 



-xlog (x) - 2(^y (|(1 -x 2 ) + x 2 log(x) 



The monomer concentration on the outside of a sphere isQ 



(X) 



„out („\ 
c sph,I\ x ) 



„out /™\ 
c sph,II\ x ) 



< x < VI + r 2 - 



C6 



; VTT 

; x > 1 



r < x < 1 



where 



c%hM = Cb(j^) (-2rx + 4 v /x (2r + x)-(2rx + x 2 ) log (l + f 



c S ph,iA x ) 



° b \4(r+x) 



'1 + 4rx + 3x 2 — 2 (2rx + x 2 ) log (x)) 



(B2) 



(B3) 



(B4) 



(B5) 



Notice that the first regimes in both cases disappear as the radius becomes large, and one 
can readily check that the second regime's profile goes to the flat case's profile in this limit. 



This concentration profile calculation, for the outside is done in ]13[ as well, where the author 
also finds these three regimes. Our result agrees with his in regimes II and (trivially) III. But 
we differ in regime I. Since his result is not continuous across the boundary between I and 
II we believe that his calculation for that case is inaccurate. Our result is continuous, and 



2 Strictly speaking this calculation assumes L <= R. 
3 Strictly speaking this calculation assumes L <= (4/3)i? 
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furthermore the full profile integrates to the full partition function for this geometry, which 
is a non-trivial consistency check. The profile for the inside is also continuous and integrates 
to the partition function for the inside of the sphere - again a non-trivial consistency check. 
The full, exact expression for the surface tension in a spherical shell is given by: 



where v(y,e) = y + e y 2 + | e 2 y 3 , and 



A 7 S, = y [ i 



v(y,e), 



(B6) 



i {8ey-Ae^/ey (2 + ey) + Ae 2 y 2 + e 2 ) + 
ey-Ae^Jey (2 + ey) + Ae 2 y 2 + e 2 f /2 



j_ 
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§y 2 + feV 
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L e 2 + |ey 3 
y - \ + ey 2 - |ey + §e 2 y 3 - |e 2 y 2 + ^e 2 



< a/1 + -2 



s/i + i-i<y< 1 

y > 1 



(B7) 



APPENDIX C: SEGMENTED ROD RESULTS 



The integral expression for the surface tension of an iV-polymer around a flat surface is: 

r N N / r 1 \ 
A 7/ = N - (1/2)" J o dz x 1] ^ [©(* " 1) + ©(I " - «<)] J (CI) 



i-l 



where = Z\ — and 6 is the step- function. This formula gives the same answers as 

j'=i 

our intuitive way of doing these integrals for (N — 1), and (iV = 2.) With sufficient patience 
or efficient programming it should be straightforward to see how the result changes as iV 
grows larger. 

The exact expression for the surface tension for (N = 2) outside of a sphere is: 



A 7 °«* = —(- 



1 + -6' 



1 \ 5 / 2 



_ (J + ^ e 2_^ e 3 + 



8 



16 



32 



(C2) 
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FIGURES 




FIG. 1. Configurational space of a rod near a surface. The rod is represented at the angle of 
contact. 
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FIG. 2. General shape of the surface tension as a function of the rod-length; L* is the maximum 
length of a line that connects any two points on the surface, without intersecting the surface. 
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FIG. 3. Rods inside a circle. 
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FIG. 4. Rods inside a sphere. 
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L sin 



FIG. 5. Rods in a flat gap - side view. 
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FIG. 6. Configurational space of a rod near a flat surface, The rod is represented at the angle 
of contact. The complete allowed angular space is obtained by rotation of the figure around the 
z-axis. 
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FIG. 7. Configurational space of rods in- and outside of a sphere. I) z > z° ut ; II) z < z° ut ; III) 
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FIG. 8. The rod concentrations, p(x)/pb, compared for rods in- and outside of a sphere and 
near a flat wall. The distances are in units of L. so: Outside of a sphere; fi: Flat; si: Inside a 
sphere. We took R = 2 L for this picture. 
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0.1 ' 0.2 ' 0.3 ~CL4 CL~5 



FIG. 9. This plot shows p out (x)/pi, weighted by the appropriate Jacobian, for the plane (solid 
line) and a sphere of radius R = L (dashed line) . The result in the text corresponds to the equality 
of the area between the solid curves to the one between the dashed curves. 
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FIG. 10. Rods in a container with an open shell. 
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FIG. 11. Segmented rods near surfaces. 




FIG. 12. This construction is used in the computation of the small L expansion of the surface 
tension. 
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TABLES 



System Spherical Cylindrical 

Quantum Particle Neutral Attraction 

Gaussian Chain Neutral Attraction 

Hard Spheres Repulsion Neutral 

Rods Repulsion Attraction 

TABLE I. Table of potentials induced by curvature in several systems. 



39 



